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In this Letter we pose the question of whether a many-body quantum system with a full set of 
conserved quantities can relax to an equilibrium state, and, if it can, what the properties of such 
state are. We confirm the relaxation hypothesis through a thorough ab initio numerical investigation 
of the dynamics of hard-core bosons on a one-dimensional lattice. Further, a natural extension of the 
Gibbs ensemble to integrable systems results in a theory that is able to predict the mean values of 
physical observables after relaxation. Finally, we show that our generalized equilibrium carries more 
memory of the initial conditions than the usual thermodynamic one. This effect may have many 
experimental consequences, some of which having already been observed in the recent experiment 
on the non-equilibrium dynamics of one-dimensional hard-core bosons in a harmonic potential [T. 
Kinoshita, T. Wenger, D. S. Weiss, Nature (London) 440, 900 (2006)]. 
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Introduction- Integrable quantum gases traditionally 
belong to the domain of mathematical physics, with little 
or no connection to experiments. However, the experi- 
mental work on confined quantum-degenerate gases has 
recently yielded faithful realizations of a number of inte- 
grable one-dimensional many-body systems, thus making 
them phenomenologically relevant. Examples include the 
gas of hard-core bosons (Girardeau) [li, i2, ^], realized in 
\^; its lattice version [5] studied in our Letter, realized 
in 6]; finite-strength s-wave- interacting spin-0 bosons 
(Licb-Liniger-McGuire) 0, Q^realized in [1, 113, El in 
the mean- field regime and in 13,ll3| in the regime of in- 



teractions of intermediate strength; and spin-^-fermions 
(Yang-Gaudin) [1J|, realized in J5j. The list has a poten- 
tial to grow to include also the fermionic p-wave version 
of hard-core particles [la uM', the gas of 1/r^ interact- 
ing atoms (Calogero-Sutherland) [ISJ, Il9l: a nd the gas of 
fermions on a lattice (Fermi-Hubbard) [20|. The experi- 
ment [6| — which is a realization of the system whose time- 
dynamics we study in the present pap er — used an optical 



lattice in the tight-binding regime [21 



. The technique 
was originally developed to reach the superfluid-Mott- 
insulator transition [23| as achieved in [24|. We should 
also mention the experimental studies [25| of a related 
nonintegrable system, the one-dimensional lattice bosons 
with finite coupling. 



other research direction concerns the effects of integrals 
of motion on the autocorrelation properties of large sys- 
tems, first studied in 271, [28( and later speciaHzed to spin 
systems 29|, [SOj ■ More recent are the studies of the onset 



of thermalization in a large quantum system [31 



in particular in a mesoscopic-size Lieb-Liniger gas [33 

The major inspiration for our work — underlying es- 
pecially Fig. 2 below — is the recent experiment on the 
non-equilibrium dynamics of one-dimensional hard-core 
bosons in a harmonic potential performed at Penn State 
University [34J . There it was found that hard-core bosons 
do not relax to the usual state of thermodynamic equi- 
librium. The question that intrigued us is whether nev- 
ertheless there exists some kind of equilibrium state 
to which a many-body integrable system relaxes in the 
course of time evolution from even a highly noncquilib- 
rium initial state — and, if so, how to predict mean values 
of physical observables in such state. 

Generalized Gibbs ensemble- We start with the lat- 
ter question, for now simply assuming that an equilib- 
rium state exists. We conjecture that then the standard 
prescription of statistical mechanics applies: one should 
maximize the many-body entropy S = k-gTr [pln(l/p)], 
subject to the constraints imposed by all the integrals of 
motion. This results in the following many-body density 
matrix: 



An integrable model possesses many nontrivial inte- 
grals of motion, and it is natural to wonder what con- 
sequences this fact may have for time dynamics and ki- 
netics. Perhaps the best known theoretical efforts in this 
vein are the attempts to explain the suppression of equili- 
bration in the Fermi-Pasta-Ulam chains by the closeness 
to various integrable models (see [26| for a review). An- 
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where {X,n} is the full set of the integrals of motion, 
Z = Tr exp[— ^^^^ A„jX,„] is the partition function, and 
{A„i} are the Lagrange multipliers, fixed by the initial 



conditions via 



fermions on a lattice: 



Tr 



irnP = (i™)(i = 0) 



(2) 



The generalized Gibbs ensemble ([1]) reduces to the usual 
grand-canonical ensemble in the case of a generic system, 
where the only integrals of motion are the total energy, 
the number of particles, and, for periodic systems, the 
total momentum. Conceptually, the ensemble ([T]) is close 
to the one E.T. Jaynes introduced in 1957 in the context 
of the so-called "subjective statistical mechanics" [35|. 
Girardeau used Jaynes's concept to study the relaxation 
of magnetization in the XY- model [36|. Below we test 
the predictive power of 111 and ([2]) on the example of 
hard-core bosons on a one-dimensional lattice, a system 
integrable via Jordan- Wigner mapping to free fermions. 
The Hamiltonian and the (quasi-) momentum distribu- 
tion of hard-core bosons on a lattice- The Hamiltonian 
for hard core bosons (HCB) on a one-dimensional lattice 
with L sites reads 



H 



T(^l^»+i 



h.c. 



i=l 



(5) 



{c,, c]} = 1, {c,;, Cj} = {c|, c]} = for all i and j. 

The corresponding fermions obey periodic (anti-periodic) 
boundary conditions for odd (even) numbers of particles: 
c^+i = (-l)*+ici. Here and below N = Y.Li^\^i = 
J2i=i KK is the particle number operator, the same for 
both fermions and bosons. 

Integrals of motion.- It is clear from the fermionic form 
([5]) of the Hamiltonian that our system possesses as many 
conserved quantities as there are lattice sites: they are 
simply the fermionic (quasi-) momentum distribution op- 
erators 



^fc = r(fc) = ^EE^-''(^)^"^'^ 



-i27rfc(i-i')/L-t 
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where 



{k, bl} = 1, (6J2 



[b], &J] ~ for all i and j ^ i; 



kt\2 



{b]y = for aU i. 



Here b^ (b]) is the annihilation (creation) operator for 
hard-core bosons, and J is the hopping constant. For our 
theoretical predictions we use a periodic lattice {bj^,i — 
bl). However, the subsequent numerical studies are per- 
formed for the more experimentally relevant hard-wall 
boundary conditions. For sufficiently large lattice sizes 
L, the difference between real physical quantities calcu- 
lated using these two settings is negligible [371] . 

Our primary observable of interest is the HCB 
(quasi-)moinentum distribution f{k) — (/(fc)), normal- 
ized to the total number of particles N, where 



fik) 
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is the HCB (quasi-)inomentum distribution operator. 

Fermi-Bose correspondence.- Our bosonic system can 
be mapped to a free fermionic (FF) one via the Jordan- 
Wigner transformation b\ 
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,17TC:, c, 



nr=ie--.'V, b, = 

Cj, where Cj {c\) is the fermionic annihila- 
tion (creation) operator. (Note that since the spatial 
density operators for fermions and bosons are equal to 
each other, c\c^ = blb^^, the inverse mapping is straight- 
forward.) Under this transformation our Hamiltonian 
([5|) becomes just the Hamiltonian for noninteracting 



where cta. (N) ^ 1 for odd N, and cta. (N) = e^'"'^^/^ 
for even N. 

Note that if expressed through the bosonic fields, the 
above integrals of motion become complicated many- 
body operators. Consider, for example, the lattice analog 
of the fourth moment of the fermionic (quasi-)momentum 
distribution I4, defined as 
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It is one of the simplest linear combinations of the inte- 
grals of motion ([5]), but it becomes a two-body operator 
in the bosonic representation. 

Fully constrained thermodynamic ensemble.- The den- 
sity matrix for the fully constrained thermodynamic en- 
semble described above reads 



PLC 



where 



Z, 



f.c. 



^f.c. '-'^P 



Tr 
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nfe(i 



e ""). The values of the Lagrange multipliers \k 
must be fixed by the requirement that the fermionic 
(quasi-) momentum distribution predicted by ([5]) be the 
same as the (quasi-)momentum distribution of fermions 
in the actual initial — or, for that matter, time-evolved — 
state of the system. This constraint leads to Afe = 

In I ~fi ' /k.\ ) ■ As we stated above, the density matrix 
given by ([5]) is assumed to predict correctly the values of 
the system's observables after a complete relaxation from 
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FIG. 1: (color online). Momentum distribution of A'^ = 30 
hard-core bosons undergoing a free expansion from an initial 
zero-temperature hard-wall box of size Ljj^ = 150 to the final 
hard- wall box of size L = 600. The initial box is situated in 
the middle of the final one. (a) Approach to equilibrium, (b) 
Equilibrium (quasi-)momentum distribution after relaxation 
in comparison with the predictions of the grand-canonical and 
of the fully constrained ((Sjl thermodynamical ensembles. The 
prediction of the fully constrained ensemble is virtually in- 
distinct from the results of the dynamical simulation; see the 
inset for a measure of the accuracy. (An animation of the 
time evolution is posted on line [4Q].) 



an initial state with the ferniionic (quasi-)nionientuni dis- 
tribution given by f^{k). Below, we test this conjecture 
numerically, using the bosonic (quasi-)nionientuni distri- 
bution as the figure of merit. 

Numerical tests- In order to verify our predictions 
we perform two series of text-book-like numerical ex- 
periments on the relaxation of an ensemble of hard-core 
bosons on a lattice from a highly nonequilibrium state. 
We have chosen to study lattices in which the final size 
L » N, i.e., the average interparticle distance is much 
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FIG. 2: (color online). Time evolution of 

the (quasi-) momentum distribution (a) and the 
(quasi-) momentum distribution after relaxation (b) of 
A^ = 30 hard-core bosons undergoing a free expansion from 
an initial zero-temperature superlattice with period four of 
half-depth A = 8J and bound by a hard-wall box of size 
L = 600, to the final flat-bottom box {A — 0) of the same 
size. The discrepancy between the result of time propagation 
and the prediction of the fully constrained ensemble (|8]) (also 
shown in (b)) is less than the width of the line. Momentum 
peaks remain well-resolved during the whole duration of 
propagation; tfi„. = 30007i/J for the subfigure (b). (An 
animation for the time evolution can also be found in [4011.) 



larger than the lattice spacing, so that our results are 
also of relevance to continuous systems [5|] • The numeri- 
cal technique has been described elsewhere J38| . 

In the first series we prepare our gas in the ground 
state of a hard-wall box, then let the gas expand to a 
larger box. For all sizes of the initial box, we find that 
the (quasi-)momentum distribution indeed converges to 
an almost time independent distribution (see Figure [T]). 
Next, we compare the result after relaxation with the 
predictions of standard statistical mechanics and of the 
fully constrained ensemble ([T]), ([5]). We find that the 
fully constrained thermodynamics stands in an excep- 
tional agreement with the results of the dynamical prop- 
agation. (See [39| for further details of the thermal algo- 
rithm.) The accuracy of the above predictions has been 
successfully verified for the whole range of available val- 



ues of the size of the initial box, from Ly 



iV = 30 



through Lj] 



600. 



In the second series (Figure [2]) we study the effects 
of the memory of the initial conditions that is stored 
in the fully constrained ensemble ^, ^. Our setting 
is very similar to an actual experiment on relaxation of 
an ensemble of hard-core bosons in a harmonic potential 
3J] . There the momentum distribution was initially split 
into two peaks. After many periods of oscillation, no 
appreciable relaxation to a single-bell distribution was 
observed. In our case the system is initially in the ground 
state of a hard-wall box with a superlattice (spatially- 
periodic background potential, see [4l| for details) with 
period 4, 
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h]b, 
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and is subsequently released to a flat-bottom hard-wall 
box T4xt = 0. Our results show that even after a 
very long propagation time, the four characteristic peaks 
in the (quasi-)momentum distribution remain well re- 
solved, although their shape is modified in the course 
of the propagation. Our interpretation of both experi- 
mental and numerical results is as follows: if the initial 
(quasi-)momentum distribution consists of several well- 
separated peaks, the memory of the initial distribution 
that is stored in the ensemble ^ prevents the peaks from 
overlapping, no matter how long the propagation time. 
Note also that the residual broadening of the peaks seen 
in Figure [21 is beyond the experimental accuracy in 34 1. 

Summary- We have demonstrated that an integrable 
many-body quantum system — one-dimensional hard-core 
bosons on a lattice — can undergo relaxation to an equi- 
librium state. The properties of this state are governed 
by the usual laws of statistical physics, properly updated 
to accommodate all the integrals of motion. We further 
show that our generalized equilibrium state carries more 
memory of the initial conditions than the usual thermo- 
dynamic one. It is in the light of that observation that 
we interpret the results of the recent experiment on the 
non-equilibrium dynamics of one-dimensional hard-core 
bosons performed at Penn State University [3J| , where an 
initial two-peaked (quasi-)momentum distribution failed 
to relax to a single-bell distribution. 
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